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5-H . 

C— . We study the phase diagram of a two component Fermi system with a weak attractive interaction. 

■^T Our analysis includes the leading order Hartree energy shifts and pairing correlations at finite 

temperature and chemical potential difference between the two fermion species. We show that 
in an asymmetric system, the Hartree shift to the single particle energies are important for the 
phase competition between normal and superfluid phase and can change the phase transition curve 
qualitatively. At large asymmetry we find that a novel but somewhat fragile superfluid state can be 
favored due to finite temperature effects. We also investigate the transition between the normal phase 

Q ■ and an inhomogeneous superfluid phase to study how gradient instabilities evolve with temperature 

Q_) ' and asymmetry. Finally, we adopt our analysis to study the density profiles of similar asymmetric 

J , Fermi systems that are being observed in cold atom experiments. 
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I. INTRODUCTION 



The superfluid nature of the ground state of symmetric two-component Fermi systems, where pairing occurs 
between equal densities of the two species, has been well established, both on theoretical and experimental 
grounds, since the pioneering work of Bardeen, Cooper and Schricffcr (BCS) fifty years ago [l|. However, the 
phase structure of asymmetric Fermi systems still remains unclear. In particular, theoretical work suggests 
several novel and competing superfluid phases may be possible. These include: (i) the Ferrel, Fulde, Larkin, 
JZZ. • Ovchinnikov (FFLO) phase, which is an inhomogeneous superfluid phase with a spatially varying superfluid 

qq | order parameter [2], y|; (ii) the breached-pair or gapless superfluid phase, which is a homogeneous superfluid 

psj . phase containing gapless fermionic excitations even at zero temperature [4|]; and (iii) phases with broken 

OQ ' rotational symmetry where pairing is facilitated by a deformation of the Fermi surface [5|, |6j. Quantum 

—J , Monte Carlo (QMC) calculations have been used to calculate the thermodynamic properties of a symmetric 

^^ ■ system at zero temperature [3] for arbitrary interaction strength. The authors of used these results to 

OO ! deduce the presence of a splitting point in the phase diagram of asymmetric systems in the strong interaction 

regime, where a gapped superfluid phase, a gapless superfluid phase and an inhomogeneous superfluid phase 
coexist. Their analysis suggests that with the interaction strength tuned to a value such that the s-wave 
scattering length between the two species is infinite, on increasing the difference in the chemical potentials 
of the two species the system undergoes a first order transition from a homogeneous gapped superfluid phase 
to an inhomogeneous phase, before eventually transiting to the normal phase. 

To investigate the possible existence of these phases, cold atom experiments are now exploring the 
thermodynamic and linear response properties of asymmetric Fermi systems with attractive interactions 
pi LUA [HI 03, [13 • These experiments which trap and cool two hyperfine states of 6 Li atoms, have unprece- 
dented control over the sample. They can: (i) magnetically tune the interaction strength between the two 
hyperfine states through Feshbach resonances; (ii) control the population asymmetry by loading different 
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numbers of atoms; and (iii) vary the temperature. In the strongly interacting regime, where the two-body 
scattering length is large, these experiments have already observed how the superfluid properties change with 
number asymmetry and temperature. Experimental measurements of the density profiles and the response to 
radio frequency probes seem to indicate that these novel phases are not realized in cold atom traps. Instead 
one finds strong indications of a first-order phase transition between a superfluid state with zero number 
asymmetry and a normal state with a large asymmetry [9|, |l2| ■ 

The absence of intervening novel superfluid phases in cold atom experiments is intriguing. These phases 
may still exist at very low temperature and in the weak coupling regime. To address if these phases can 
be realized in experiments in the future, we investigate the phase structure of asymmetric Fermi systems 
at finite temperature and establish the parameter region where interesting new phases of superfluidity can 
be realized. Our study differs from similar investigations reported in Ref. [14 , Il5l | in two ways. (For a 
discussion of the phase structure of asymmetric systems at finite temperature in condensed matter systems, 
see [IflllsllMllyl-) Fhst, we restrict our analysis to the weak coupling region and perform a self-consistent 
calculation of the thermodynamic properties of both the normal and the superfluid phases within the purview 
of mean field theory (Hartree approximation) . One drawback of this approach is that we neglect particle- hole 
screening which is important in the gap equation at weak coupling [2fJ, [21| . This screening is expected to 
reduce the gap by a factor ~ 2.2. In the conclusions (section IIV[) we will return to a discussion of how 
this suppression may affect our final answers. Second, we establish the region in temperature and number 
asymmetry where the homogeneous states are unstable with respect to small amplitude perturbations. 

Although our analysis is strictly valid only in the weak coupling regime, unlike earlier work, we properly 
account for the single particle energy shifts - or Hartree corrections - that are present both in the normal 
and superfluid state, in a self-consistent field theoretic approach. We find that these corrections can change 
the location and shape of the phase boundaries on the phase diagram of the two component Fermi gas. 
We also find that a new gapless, finite temperature superfluid state is favored over the normal state at 
temperatures above the critical temperature for the first-order transition. In a narrow window of chemical 
potential asymmetry, with increasing temperature there is first a first-order transition from the superfluid 
to the normal phase and then at higher temperature a second-order transition to a weakly superfluid state 
which exists in limited temperature interval. In this state, finite temperature effects that smear the Fermi 
surface, facilitate pairing. We will refer to this homogeneous phase as the "fragile" superfluid and it has 
lower energy than the normal state but is unstable with respect to gradient instabilities. 

To study these instabilities we expand the free energy in a Ginzburg-Landau series about the solution of 
the gap equation, and calculate the coefficient of the quadratic term in the order parameter as a function of 
its Fourier mode index. We ask the question whether this coefficient is negative for some Fourier mode(s), 
indicating a gradient instability. The simplest case is an expansion about the normal phase with a zero 
background value for the difcrmion condensate. In the presence of a small position dependent condensate 



, ? (r) = / S 3 ke lkr ?? (k) , (1) 

we write the change in the free energy as, 

J a 3 k(a + /(|k|))r ? (k)r 7 (-k) + Otf) , (2) 

where we have split the terms in a way that /(|k|) is zero for k = 0. a < points to an instability 
towards the formation of a homogeneous condensate while /(|k|) < points to an instability towards an 
inhomogeneous modulation of the condensate. We find that there is a window of temperatures and chemical 
potential splitting where the inhomogeneous superfluid phases may be favored. 

The plan of the paper is as follows. We begin by writing down the lagrangian, and reviewing how to 
include Hartree corrections to the calculation of the gap parameter A and the free energy Q of the system 
in section [Hi We then proceed with the calculation of A and il in section Hill which will let us find the phase 
boundary between the superfluid and the normal phase. We first consider homogeneous phases, meaning 
A independent of position, and follow with a discussion of inhomogeneous phases. We conclude with a 



summary of our results in section HVl where we use our expressions to calculate the polarization as a function 
of the distance from the center of an isotropic atomic trap, for typical trap parameters. We also identify the 
region where inhomogencous phases may be found in a typical trap geometry. 

II. THE MODEL LAGRANGIAN 

We describe the gas of two species of fermions, ipi and ip 2 , with chemical potentials /ii and fj, 2 , at finite 
temperature T, by a model lagrangian density of the form, 

C = ii [iidt - £{p))5 af 3 + Sfialpj i/>p + -i>\i\ijii a , (3) 

where £(p) = p 2 /(2m) — fi. The chemical potentials for the two species of fermions, i>\ and 1^2, in terms of 
the average chemical potential fi and the splitting 26 fi arc, fi\ = /i + 5fi and [i 2 = ^ — Sfi, respectively. We 
are interested in the case where the interaction between the two species of fermions is attractive, meaning 
A > 0. In the BCS regime, then, for small enough T and (5/x, the phase of the system will be characterized 
by a non-zero difermion condensate 

(Mr)Mr)) = ^A(r) . (4) 

We are interested in the phase boundary between such a superfluid phase and the normal phase, where 
A(r) = 0. 

In addition, the system will be characterized by specific profiles for the number density of the two species , 

(4 (rtyi (r)> = m (r) , (4 2 (r)0 2 (r)} = n 2 (r) . (5) 

In our calculations we will find it more convenient to write n\ and n 2 in terms of the average density 
n(r) = (l/2)(ni(r) + n 2 (T)) and the difference in densities Sn(r) = (l/2)(ni(r) — rz.2(r)). 

To specify the system for any T, /i, Sfj, for a given A (which we will trade for the scattering length a), we 
need three equations that let us solve for the three variables A, n and Sn. These are the gap and the number 
equations. For homogeneous condensates these take the form, 

in 

-26n . (6) 

We calculate the free energy, fi, in a mean field approximation where we replace the four Fermi interaction 
by its mean field value 
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+Ani(r)^|(r)^ 2 (r) + \n 2 (r)ip\(r)ii(r) - Arii(r)n 2 (r) . 

The terms proportional to ri\ and n 2 give rise to Hartree corrections to the free energy 

Upon making the mean field approximation, and performing standard manipulations, we can write the 
lagrangian density (eq. ([3])) in a quadratic form in terms of the Nambu-Gorkov spinor, 

* = (Vi i2 4 4 2 ) T ■ (8) 

The final answer is 



where £(p) = p 2 /(2m) — ji. jx and <5/2 include the Hartrec terms and are given by 

ji = n + Xn , 6 p. = <5/i — XSn . (10) 

That the Hartree corrections to the chemical potentials should be this way can be seen immediately if we 
recall that the Hartree terms change \x\ to ji\ = fii + \ri2 and Li2 to p,2 = M2 + Ani. The presence of the 
mysterious looking term S(0)jl can be explained as follows. 

When written in terms of the fields ip and ip\ the mean field lagrangian density has a piece ^ Q JAaipJ^ipa. 
To write this in a symmetric form in terms of the components of the Nambu-Gorkov spinor 'J, we need to 
exchange the ordering of i/j and ip\ which gives rise to the term in question. More explicitly 

Y^j MaV'aV'c = 2 Yi^alpi^a + Aa^a^a) = £ ^ifia'^a ~~ AaV'aVi + fl a S(0)) , (11) 

a a a 

where we have used the fermion anticommutation relation 

{^ a (r),^(r')} = 6^8(r-r'). (12) 

This term is important in canceling out a divergent contribution to the free energy, as we shall see below. 
A similar term occurs while reordering ^2 a tj)' a {idt — p 2 / (2m))ip a , but is the same in normal and superfluid 
matter and does not affect the phase competition. 

In the mean field approximation the lagrangian density is bilinear in the fermion fields and the free energy 
is found by direct integration over the fields. We find 

I d 4 x E n = fd 4 x E {-S(0)fi + ^^p + A(n 2 (r) - 6n 2 {r))} 

-{[5**(-**:i$ + " s^ + s,)^-w}- <13) 

where xe represents the euclidean space four vector (a; 4 ,r), with r lying in position space with volume V 
and a; 4 £ [-1/(2T), 1/(2T)]. 

We now compare the free energies of the superfluid phases to the normal phase to find where the boundary 
between the phases lies. 

III. PHASE BOUNDARY BETWEEN THE SUPERFLUID AND NORMAL PHASE 

We first consider the competition between the homogeneous superfluid phase and the normal phase, and 
then the competition between inhomogeneous supcrfluids and the normal phase. 

A. Homogeneous superfluid and the normal phase 

The homogeneous phases are defined by the condition that A, n and Sn are all independent of r. The 
argument of the log in eq. (|13p is then diagonal in momentum space and the free energy density is simply, 



n(A,n,Sn,T,n,6Li)= {-ju /# 3 p(l) + J-J- + A(n 2 - 6n 2 )} 

-{/ S3 P[| E log[(tp 4 + Sp-e(p))(2 P 4 + Sfi + i(p)))+(dp^-dp)]} 

p 4 = 
{2n+l)TvT 

-/ 3 3p{T[lo g [co S h(^)] + log[cosh(^±M)]] +A } , 



(14) 



where e(p) = V(£ 2 (p) + |A| 2 ) and we have rewritten 6(0) as an intcrgal over momentum space. 

The gap equation is the condition that the free energy is stationary with respect to small variations in 
the magnitude of A. (A position independent phase of A does not affect the free energy and hence, for 
simplicity, we will take A to be real and positive for the rest of the section (|III A[) .) 
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The trivial solution, A = 0, corresponds to the normal phase. A competing supcrfluid phase exists if eq. (|15p 
has a solution with A^O. 

As it stands, the second term in eq. (|15p is linearly divergent in p and needs to be regulated in some 
manner. One way is to cut off the momentum integration at some momentum A, chosen to be sufficiently 
larger than the Fermi momentum, kp — y/(2m[i), to capture all the features of the intergand. In the weak 
coupling (BCS) regime, it enough to take A to be several times kp. The solution of eq. (fT5)) then depends 
on A and A, but it is useful to rewrite the results in terms of a physical observable. A popular choice is to 
use the relation between the s-wave scattering length a, and A. 
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Canceling 2A ^ from the right hand side of the gap equation, eq. (fT5"|) . and substituting 1/A in the first 
term in eq. (115p from eq. (|16[) . we obtain the relation 
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which is now cutoff independent. Recall, however, that A also appears implicitly in the definitions of p, and 
5ji (eq. (fTU|) ) where it multiplies n and 6n, respectively. In all these places, we replace A by — (Aira)/m. To 
see why this is reasonable, consider eq. (|16l) with a momentum cut off A. 
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If |47rAa| -C 1, then in eq. (p~8|) we can ignore the second term on the right hand side compared to the first 
term. In the weak coupling limit |7rfci?a| « 1. If A is taken to be not many times larger than kp, as we 
argued can be done in the weak coupling regime, then indeed we can take A to be — (Ai:a)/m to a good 
approximation. 

To summarize, the gap equation is given by eq. (|17p with 
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(19) 



The number equations are found by explicitly calculating the derivatives with respect to fi and 6fi (eq. ©). 
The final expressions are given below. 
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(21) 



with fi and Sfi given by eq. (|19p . 

For the normal phase, A = 0, and we can solve eqs. (|2QI2ip to find n and Sn to obtain values we will call 
njv and Sn^ respectively. The free energy of the normal phase is given by, 



Q 



N (n,6n,T) = j( )(n 2 N -Sn 2 N )j 



(22) 
-/s3p{T[log[cosh(^(Hk) ] + log[cosh(^^±i^)]] + » N } 

with £(p)jv = p 2 /(2m) — jx^ and 

M« = M + ( ) n N , Sp, N = Sn - ( ) Sn N . (23) 

V m I \ m ) 

If, in addition, eqs. (|17I20I21[) possess solutions with A ^ 0, ua and Sua, we need to compare the free 
energies of these supcrfluid solutions to iljy. The difference in the free energies is given by 

-/^p{T[log[cosh(^±i^)] + log[cosh(^±!(P))]] +,-^} . 

(24) 

If Q s — flpf > then the normal phase is favored over the supcrfluid phase, and vice versa. In Fig. |T]) we 
look at the boundary marking the normal to supcrfluid phase transition in T, 8[i space for g = (irkpap -1 ' = 
—0.72. For different values of u and m, if we choose a so that the dimcnsionlcss parameter g remains the 
same, then the physical quantities scale with [i and m as follows. 

n a {m, m, T, <5 M ) = m^V 3 / 2 >n Q (l, 1, (T/jk), (S/i/n)) 
A(m, n, T, 6/i) = A*A(1, 1, (T/ M ), («*/////)) (25) 

n(m, M , r, Sn) = m( 3 / 2 V (5/2) fi(i, 1, (t/a*), (<Wm)) • 

In particular, the value of Ao, the value of the gap parameter at T = and 5/x = 0j scales linearly with \i. 
To remove the dependence of the phase transition curve on the overall scales, it is useful to draw the phase 
diagram in terms of dimcnsionlcss variables. Since Ao is proportional to /i, we take the x — axis to be S/j,/A.q 
and the y — axis to be T/Aq. 

For comparison, in Fig. §T§ we also show the result when one does not include the Hartree corrections. One 
effect of Hartree corrections is to simply shift the chemical potentials. For example, at T = and Sfi = 0, 
including these increases the "effective" chemical potential, jj, = /j, — (Aira/m)n. (Recall that a < in the 
BCS regime.) Therefore, for the same a, the value of Ao is greater when we include the Hartree corrections, 
compared to when we do not include them. To get rid of this overall change in Ao , we scale T and Sfi in 
the "non-Hartree" curve by the "non-Hartree" value of Ao. The two curves are still different, and that has 
to do with the fact that including Hartree corrections affects the competition between the normal and the 
supcrfluid phases. Number densities, and therefore Hartree corrections, are different in the two phases at the 
phase boundary, if the phase transition between the two phases is first order. In the following paragraphs, we 
discuss the effect of Hartree corrections on the phase transition curve between the normal and the supcrfluid 
phase in more detail. 

First looking at T = 0, S/j,/Aq at the first order transition is larger than the weak coupling value 0.707, as 
pointed out earlier in [7j. Since at T = 0, Sn = in the superfluid phase (the supcrfluid phase is gapped), 
this is not simply due to a reduction in the effective splitting between the Fermi surfaces in the superfluid 
region (Sfi = Sfx ~ O.8IA0 for the superfluid phase at the phase transition). As discussed above, the change 
is due to the fact that n and Sn change abruptly at the first order phase boundary. More specifically, there 



arc two effects both of which drive the phase transition to larger Sfi. Firstly, Sn is positive in the normal 
phase, implying Sfl is smaller than 5fi, which increases the free energy of the normal phase. Secondly, n is 
larger in the supcrfluid phase because of pairing, and this also makes the supcrfluid phase more favorable. 
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FIG. 1: (colored online) We show the phase boundary between normal and homogeneous superfluid phases in 5/j,, T 
space, for g = —0.72. In the curve marked by circles (black online), we include the Hartree corrections. The value of 
Ao/fi is 0.058. The curve marked by squares (red online) does not include Hartree corrections. The value of Ao/fi 
in this case is 0.03. At T = 0, the first order phase transition from the superfluid to the normal phase occurs at 
Sfi/ Ao ~ 0.71, consistent with weak coupling results. The filled points mark the region where the phase transition 
is first order while the hollow points mark a second order transition. The points where the transition changes from 
being second order to becoming first order have been encircled. 



Next, looking along the y — axis, at 6fi = 0, the transition from the superfluid to the normal phase at the 
critical temperature is second order, and the number densities in the two phases arc the same at the critical 
temperature. It is understandable, therefore, that the transition occurs at T c /Aq ~ 0.567, the standard 
weak coupling value. 



Things are interesting close to 5[i/Aq ~ 0.8, where the shape of the curve is qualitatively altered. For a 
window of splittings, <5/i/Ao € (0.79, 0.81), as we increase the temperature, we encounter not one but three 
normal-supcrfluid transitions. To clarify how this comes about, consider the shape of fi as a function of A 
for various values of the temperature, at 5fi/Ao = 0.795 (Fig. [5]). At T = 0, the local minimum at A = Ao 
is favored. As we increase the temperature, this minimum becomes shallower and eventually there is a first 
order transition to the normal phase. As we keep increasing the temperature, the A = solution becomes 
unstable and the superfluid state is favored again for a range of temperatures. Eventually, at a large enough 
temperature, there is a second order transition to the normal phase. This reappearance of superfluidity at 
higher temperatures can be understood intuitively as follows. At zero temperature, BCS pairing is stressed 
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FIG. 2: (colored online) This figure shows the free energy as a function of A at S/i/Aq = 0.795 for various values of 
T. The solid curve (black online) corresponds to the lowest temperature, T/Aq — 0.003, and clearly shows a robust 
superfluid phase with A ~ Ao. The dashed curve (blue online) corresponds to T/Ao = 0.08 and is the point of a 
first order phase transition to the normal phase. A "fragile" superfluid phase reappears before it disappears again at 
a higher temperature. 



due to a non-zero S/i, because fermions of the two species can not find partners of opposite momenta lying on 
the distinct Fermi surfaces determined by their different chemical potentials. BCS pairing in such systems 
requires the two distinct Fermi spheres to equalize at a Fermi momentum different from the value given by 
the corresponding chemical potentials, costing free energy. At non-zero temperatures, however, the Fermi 
surfaces are smeared and it is possible to find partners of opposite momenta even without equalizing the 
Fermi surfaces. Such an effect of temperature on p airing has been noted previously in the context of pairing 
between u and d quarks in the 2SC phase [22, 123, 0, 12a, 123, 1221, l28( , and a nuclear physics context [29| ■ 

To see in a different way why there are three transitions, it is useful to expand the free energy in powers 
of A 2 (Ginzburg-Landau expansion) and look at how the coefficient of the quadratic term changes as we 
change T, keeping Sfj, constant. We write the free energy as aA 2 + C(A 4 ). a < points to an instability 
in the A = state to the formation of a non-zero condensate, while a > means that the normal state is 
locally stable, but does not tell us whether it is globally favored or not. 

In Fig. ([3]) we plot a (in units of A 2 ,) as a function of T/Aq, for <5/i = 0.795Ao kept constant. For small 
T, the normal phase is locally stable, but globally disfavored to the A = A state (Fig. @). As we increase 
T, the superfluid phase becomes less and less favorable until finally at T/Aq = 0.08 the normal phase is 
globally favored. On increasing T further, we find that the normal phase becomes locally unstable in the 
region T/Aq € (0.22,0.33), and this gives rise to a region of "fragile" superfluidity. 



For smaller couplings the window of S)i/Aq where we encounter three normal-supcrfluid transitions, be- 
comes narrower, while for stronger couplings it becomes wider. 
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FIG. 3: Plot of a as a function of T/A for 8fi = 0.795A . The normal phase is locally stable for T/A £ (0,0.22), 
locally unstable for T/Ao £ (0.22, 0.33), and locally stable for T larger than 0.33Aq. 



B. Inhomogeneous condensates 



Let us now consider the case where A, n and Sn depend on r. The argument of the log in eq. (|13p is no 
longer diagonal in momentum space and hence we evaluate the free energy in a Ginzburg-Landau expansion. 
We are interested in finding the curve along which the normal phase becomes unstable to the growth of an 
inhomogeneous condensate. Working in the limit of small A(r), we expand the log in A(r) and drop terms 
proportional to A 4 and higher, and obtain 



f d 3 ril= [ d 3 r{-jl f a 3 p{l) 



|A(r) 



X(n 2 -Sn 2 )\ 



{T J2 [Tr r log(^ 4 + ^-e( P ))(*P 4 + ^ + |(p))) 



p 4 = 

(2n + l)7rT 



(26) 



- Tr r ((ip 4 + 5p - i(p))^A(r)(ip 4 + Sfi + l(p))^ A* (r))] } 



We now argue that within the approximations we are working in, we can replace n and Sn by their values 
in normal matter to compute Q. It is obvious that we can do so in the term, 



T J2 ((*P 4 + ^-e(p)) ( - 1) A(r)( l p 4 + ( 5 / i + e(p)) ( - 1) A*(r)), 



(27) 



p 4 = 

(2n + l)7rT 



because the corrections to n and Sn due to pairing is proportional to A 2 , and keeping these corrections 
in eq. (|27|) (where they appear in p, and Sfi) will only change the result by order A 4 . There are C(A 2 ) 
contributions to fi from /2/S 3 p(l), X(n 2 — Sn 2 ) and Tr r log((ip 4 + 5 p. — £,(p))(ip 4 + Sp + ^(p))), but in all 
these cases, the A 2 correction to n or Sn is further multiplied by A and for weak coupling, this should give 
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a small overall contribution. (This is the A 2 correction to the Hartree term which itself is a correction) . In 
the approximation where we neglect this contribution, when we calculate the difference between the normal 
and supcrfluid free energies, these three terms cancel out and we obtain 

jd 3 r(n(A(r)) - Q N ) = J d 3 r{ 1 -^-} 

+ h 5 v{T Yl [^r((i/ + ^iv-|(p)iv) ( - 1) A(r)(ip 4 + ^iv + |(p)jv) ( - 1) A*(r))]} 

(2ti. + 1)ttT 

(28) 
In momentum space this gives, 

[ d 3 r(n(A(r)) - n N ) = f 5 3 kA(k)A*(-k)[-^+ 

+ /> p {™+t ]T ^ 4 + ^ Ar -|(p + k) J v)(- 1 )^ 4 + 5 / i JV + f(pV) ( - 1) }] , 

J Mr 4 

(2ti+1)tvT 

(29) 

where we have used eq. (|16[) to rewrite A. 

It is convenient at this point, to separate the "potential" contribution (the contribution independent of k) 
from the "gradient" contribution (zero for k = 0) [3fj|. We write 



with, 



J d 3 r(n(A(r)) - n N ) = J 3 3 kA(k)A*(-k){a + f(\k\)} 

-£L + y> p {™+ T J2 (*/ + ^-f(pM ( - 1 W + ^ + £(pM ( - 1) } 



(30) 



(2n + l)7rT 



m / _.n r m 



tanh 



■6m + £(p)N\ . , ,/-tim + i(p)N- 



(31) 



2T 



tanh(- 



2T 



4£(p) 



and. 



/(|k|) = i />p{r £ 

i 



K( P +kv-e(p) 



iVj 



((zp 4 + 5~ m y - (^(p)at) 2 )(^ 4 + <^a0 2 - (£(P + k) w ) 2 ) ■ 

(e( P + k)^ - e(p)iv) rsdCpV) <?(I(p + k)*) 



(2~+l)iT 



(32) 



(ftp + k^+^p)*) 1 C(P) 



AT 



where 



fl (£) = I[tanh(fc^)+tanh v 



e(p + k)jv 

-dfi + i, 



}■ 



(33) 



At any given temperature, for large enough 5/x, the normal phase will be favored over a phase with non- 
zero A(k), and the combination a + f(k) will be positive for all values of k = |k|. As we decrease Sfi keeping 
T constant, a + f(k) may become zero, and then negative, for a single mode with momentum k = fc m i n . 
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FIG. 4: (colored online) We show the curve marking a second order phase transition between normal and inhomoge- 
neous superfluid phases in 8fi, T space, for g = —0.72. In the curve marked by circles (black online), we include the 
Hartree corrections while the curve marked by squares (red online) does not include these. Where the symbols are 
hollow, the value of k at which a + f(k) becomes zero, is zero, meaning that the instability is towards the formation of 
homogeneous condensates. It is reasonable that at smaller Sfi and larger temperatures, the formation of homogeneous 
condensates is favored because the Fermi surface is smeared out and it is no longer advantageous to form Cooper 
pairs with non-zero net momenta. 

If fc m i n 7^ 0, this point symbolizes the onset of the instability towards the formation of an inhomogeneous 
condensate. At lower S/j,, more momentum modes may become unstable. If the transition from normal to 
inhomogeneous superfluidity is actually first order, then we expect it to occur for values of 8fi larger than 
the value we find using this second order analysis. 

Fig. Q shows the curve in T, S/j, space, which tells us the value of S/j, where the coefficient a + f(k) 
becomes zero for some k, as we decrease Sfx from a large value keeping T constant. Since we are looking at a 
second order phase transition line, once we take into account the shift in the chemical potentials due to the 
Hartree corrections, the results are consistent with the well known results for weak coupling. For example, 
at T = 0, the value of Sp, = 6/J, — (—4:ira/m)8n at the phase transition is given by 5j1/Aq ~ 0.75. Weak 
coupling Ginzburg-Landau calculations tell us that the transition should occur at Sfi* = 0.754Ao with the 
most unstable momentum fc m ; n given by Wffc m i n /2 = 1.2<5/i*. We find that the value of k which becomes 
unstable at <5/2/A ~ 0.75 satisfies vpk m i n /2 ~ 1.245/2. For reference we also show the LOFF boundary 
excluding the Hartree corrections using square dots (red online). The instability towards inhomogeneous 
condensates at T — develops at Sfj,/ Aq ~ 0.75, and the value of vpk m i n /2 we find at this point is 1.37S[i. 



IV. CONCLUSIONS 



We calculate the effect of Hartree corrections on the phase transition curve separating the normal from the 
superfluid phase in the T, 5fi plane, in the BCS regime. Our analysis is similar in spirit to the calculations 
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FIG. 5: (colored online) We show the phase boundaries between the normal and the inhomogeneous superfluid phases 
and the normal and the homogeneous superfluid phase, for g = —0.72. The lower curve (red online) represents the 
homogeneous superfluid boundary, while the upper curve (black online) marks the instability towards inhomogeneity. 
The two curves coincide for S/j, < 0.77Ao. The shaded region is where LOFF-like phases may exist. The three parallel 
thick lines (blue online) placed at 5/i/Ao = 0.750, 0.770 and 0.795, show the stability of the mean field background 
solution to the growth of inhomogeneities. In the temperature range for which a line is solid, the homogeneous 
condensate is stable, while for the temperatures for which a line is dotted, the condensate is unstable. The input for 
these lines comes from the data shown in Fig. (0. 



performed by the authors of [Ij|, who however concentrate on the unitary regime. In their work, the Hartree 
corrections are taken to be of a well motivated form [3l[ chosen so as to reproduce the data from Monte 
Carlo simulations performed at unitarity Q- We work in the BCS regime, where the Hartree corrections 
can be calculated from first principles, as described above. Furthermore, wc consider the instability of the 
homogeneous states with respect to the formation of inhomogeneous condensates. 

The results are shown in Figs. (JTJ [3} for g — —0.72. In figure Fig. ©, we zoom in on the large S/i region 
and show the homogeneous-normal and inhomogcneous-normal phase boundaries on the same diagram. The 
curve that is non-zero up to Sfi ~ 0.99Ao (black online) shows the parameter values for which the coefficient 
of the A(k) 2 term, in a Ginzburg-Landau expansion about the normal phase, becomes negative for some 
value of k = fc m i„. This signals an instability towards the formation of an inhomogeneous condensate. The 
curve that meets the x— axis at 8fi ~ O.8IA0 (red online) represents the points where the free energies of 
the homogeneous superfluid and the normal phase are equal. Therefore, the region between the boundaries, 
highlighted by a shading (gray online), is the region where we expect LOFF-like phases to exist. Our first 
conclusion is that at weak coupling, Hartree corrections do not destroy the parameter space where LOFF-like 
phases may exist. 

Secondly, we find that the shape of the curve marking the transition from the homogeneous superfluid 
phase to the normal phase, is qualitatively altered (Fig. {1])). For a narrow range of 8[i near the first 
order transition, we encounter three normal-supcrfluid transitions as we increase the temperature keeping 
Sfi constant. There is a first order phase transition from the superfluid to the normal phase at a low 
temperature. As we increase the temperature, a "fragile" superfluid phase, where pairing is assisted by a 
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temperature induced smearing of the Fermi surfaces, reappears in a window of temperatures. 

As we shall see below, analysis of fluctuations about the mean field will tell us that this "fragile" super- 
fluid may give way to an inhomogeneous phase (Figs. ([5l [7])). For 8/i/Aq = 0.750, 0.770 and 0.795, the 
temperatures for which the homogeneous condensate is unstable is shown by the dotted lines (blue online) 
in Fig. ©. 

If this mean field picture is taken as it is (and ignoring the inhomogeneous phases for a moment), this 
structure of the phase diagram could be observed in cold atomic traps. As we go away from the center of 
the trap, the effective chemical potential and therefore Ao decreases. Since T and S/i are constant across the 
trap, this implies that T/Aq and <5/i/Ao increases. The "fragile" supcrfluid region can be probed by tuning 
the parameters, namely the number of the two species in the trap N± and N 2 , the scattering length a, and 
the temperature T, so that for a given trap geometry, we pass through the "fragile" superfluid region as we 
go from the center of the trap to the outside. 




FIG. 6: (colored online) Polarization as a function of the r plotted in units of n. The bold line (black online) 
corresponds to T = 0.15 fi and the dashed line (red online) to T = 0.365 /i. Also shown are regions where the system is 
unstable to the growth of inhomogeneities. The thin continuous lines (black online) specify this region for T — O.lS/i 
and the thin dashed lines (red online) mark this region for T — 0.36S/J,. The lower boundary marked by a thin dotted 
line, coincides with the left boundary of the thin continuous line. The upper boundary, however, appears at a larger 
r for the continuous line, telling us that the window of inhomogeneous pairing is wider at lower temperature. For 
T = 0.16/1 the total value of N in the trap is 3.21 x 10 (+7) and the value of 6N is 1.26 x 10 (+6) . For T = 0.365/i the 
total value of N = 3.20 x 10 (+7) and the value of SN = 1.53 x 10 (+6) . 



In Fig. ©, we plot profiles of the polarization p = (n\ — n^/ini + n 2 ) as a function of r, where r is the 
distance from the center of the trap and ri\ and n 2 are the number densities of the two species. The t rap 
is taken to be spherically symmetric, with a harmonic trapping angular frequency ui = 100 rads/s [9j, llOi . 
[ill fl2l . fl3| . The harmonic potential can be written as V(r) = 1/2 x Tiu) x (r/ro) 2 with tq = y/(h/(mu>)), 
where m is the mass of 6 Li atoms. In natural units, u> = 2.599 x 10' -10 - 1 (eV), m = 5.61 x 10' +9 - ) (eV) and 
r = 0.82816 (eV)" 1 . 
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With the geometry of the trap specified, we now try to tune the parameters of the experiment N = 
(Ni + N 2 )/2, 5N = (Ni — N 2 )/2, a and T and try to identify the region where there is a phase transition 
from the superfiuid phase to the unpaired phase. This is interesting because of two reasons. Firstly, since 
the polarization rises rapidly at this phase transition boundary, this boundary can be recognized in the trap 
by observing the polarization as a function of r. Secondly, the inhomogeneous phases are likely to be found 
near this boundary at low temperatures, and we would like to estimate the range of r where these phases 
are likely to exist. We choose a = —0.0205 (eV)" 1 , for which the atomic system, especially at the center of 
the trap, can not be considered weakly coupled, and hence our approximations may not be quantitatively 
accurate near the center of the trap. But the large coupling amplifies the features in the polarization profile 
that we are looking for. We use these features to deduce qualitative conclusions which may be more general 
than our approximations. 

Although the true experimental variables are N and 5N , we find it more convenient to fix the values of the 
average chemical potential at the center of the trap, fi$, and the chemical potential splitting <5/i, and calculate 
N and SN in terms of these. The average chemical potential at a point r in the trap is jx{r) = /.to — V(r). 
Plotted in Fig. © is p as a function of r for fi = 1.314 x 10(~ 7 ) (cV) and Sfi = 1.2 x 10( _8 ) (cV) for two 
different temperatures, T = O.lSfi and T = 0.36(5/i. We focus on the region near the place where we observe 
the jump in the value of the polarization. The radius r has been plotted in units of r\ defined as the radius 
where the effective chemical potential of i\)\ species is zero, i.e. r\ = ro^/(2(/^o + 5/j,)/oj) = 33.22ro. 

We see that, as expected, the jump in the polarization at the phase transition boundary becomes less sharp 
as we increase the temperature. We also mark the region where the system is unstable to the formation of 
inhomogeneous phases. We also see that the range of r where inhomogeneity may develop, decreases in size 
as we increase the temperature. 

But we expect that this mean field picture will be modified due to quantum corrections. One thing to 
note from the doublc-dot-dashcd curve (red online) in Fig. (|2|) is that the free energy curve as a function of 
A is very shallow. The mass of the fluctuations in the magnitude of the condensate is therefore very small 
and therefore we expect quantum corrections to be important in this "fragile" superfiuid state. 

The role of quantum corrections can be understood by considering the effective action describing the 
fluctuations of the two fermion condensate about the mean field value. If the system is stable to small 
fluctuations of the condensate, we may calculate their contribution to the free energy and see if they increase 
the mean field value of the free energy or decrease it. Recently, the effect of fluctuations on the free energy 
function has been calculated for T = and 5fi — [32]. These calculations suggest that the mean field 
calculations overestimate the value of Aq as well as the value of n in the superfiuid state. 



The first step in such a calculation is to find the effective potential describing the dynamics of the fluctu- 
ation field and seeing if the system is stable to small fluctuations in the mean field value. In the following, 
we analyze the stability of the "fragile" superfiuid to fluctuations that depend on position and are purely 
imaginary for real A. I.e., we consider a condensate of form 

(^ a (r)^(r)) = -e a0 (A + irj(r)) , (34) 

where A and r\ are both real. The reason for considering purely imaginary fluctuations (eq. (|3~4"|) ) is that it 
is clear from the free energy curve at T = 0.24Ao in Fig. |(3J) that the free energy is at a local minimum at 
the mean field value A, meaning that the system is stable to spatially uniform variation in the magnitude 
of the condensate. One could potentially consider position dependent fluctuations in the magnitude of the 
condensate but we only consider real rj(r) in eq. (|34p . which is related to position dependent fluctuations 
in the phase of the condensate. The effective potential describing the dynamics of r\ can be found in a 
Ginzburg-Landau expansion in the same way as was done in section Till Bl where we considered inhomoge- 
neous fluctuations about the normal phase. The only difference now is that the background has a non-zero 
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FIG. 7: /min as a function of T at three values of Sfi, Sfi = 0.750A , Sfi = 0.770A and 8(i - 0.795A () . We take 
g — —0.72 as before. Looking first at S/j, = 0.795Ao, / m in is clearly negative in the unpaired region, corresponding 
to T/Ao £ (0.08,0.22). It is also negative in the "fragile" superfluid region lying in T/A £ (0.22,0.33). Now 
considering Sfi = 0.770Ao, we find that the homogeneous condensate is unstable in T/Aq £ (0.16,0.29). Finally for 
Sjj, = 0.750Ao (and lower), we find that the homogeneous condensate is stable for all T. These three curves form the 
basis of the thick dashed lines (blue online) in Fig. ((5]l . 



condensate value, A. The final answer has the form, 



d 3 r(n(A(r)) - n N ) = y a 3 k ?7 (k)77(-k){/ A (|k|)} 



ith, 



/A(|k|) = - 



(t 



1 / I? 1 

P 17I( 



(f(p + k) - |(p)) 2 r 5 (c(p)) g(e(p + k)) 



e(p + k) 2 -e(p) 2 )L e(p) e(p + k) 



}■ 



(35) 



(36) 



The form in eq. (|35[) closely resembles the form in eq. (|30[) . except now a — 0, which can be easily understood 
as the consequence of the fact that a position independent fluctuation in the phase of the condensate can 
not change the free energy. 

In Fig. (|7|) we plot the minimum value of /a(&) (minimized over k = |k|) as a function of T at 5/i/Ao = 
0.795, 0.770 and 0.750 for g = -0.72. For <5^/A = 0.795, we see that the minimum value, f{k min ), is 
negative in the "fragile" superfluid region lying in T/Aq £ (0.22,0.33), indicating that it is unstable with 
respect to developing a position dependent phase modulation. We see (Fig. 0) that as we decrease 5fi, the 
homogeneous superfluid region becomes stable for all T, as it should. 

In light of this instability of the "fragile" superfluid, our results of the polarization profiles for the atomic 
trap may need to be modified. To find the correct profile, we will need to calculate the number densities 
of the two species in an inhomogeneous condensate more accurately than we have done above. We leave 
this calculation for future work. Another issue that warrants further investigation is the role of polarization 
effects (also known as the Gorkov-Barkhudarov corrections) that can reduce the gap. In weak coupling, 
medium polarization results in screening the interactions and reduces the gap and the critical temperature 
by a factor of (4e)^/ 3 ^ ~ 2.2 when S/j. = [20,|21(. If we were to incorporate this reduction by modifying the 
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four-fermion coupling in the pairing channel to mimic the Gorkov suppression at T — and <5/i = 0, it is easy 
to argue that the phase boundary (shown in Fig. [1]) would simply be scaled by the same reduction factor 
2.2 both in temperature and Sfi. However, a systematic calculation of these medium effects at moderate 
coupling is challenging and preliminary investigations suggest that the gap is not as strongly suppressed 
[33| |. Further, even in weak coupling the nature of this suppression is not understood at finite population 
imbalance. We leave a systematic calculation of how Gorkov corrections evolve with coupling strength and 
Sfi for future work. 

Despite these aforementioned caveats, we reiterate that our investigation here points to several interesting 
qualitative features that arise only when Hartree corrections are included. In particular we see that the 
region in phase diagram susceptible to gradient instabilities is moderately enhanced by these corrections. 
The shape of the normal-superfluid phase boundary is shown to depend on the nature of these mean field 
energy shifts since they differ in the two phases in the vicinity of the first-order transition. In particular 
the Clogston-Chandrashekar point, which is given by Sfi = A/\/2, is shifted to higher values of Sfi by these 
corrections. 
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